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Abstract 

We investigate the spectrum and eigenstates of ultracold fermionic atoms in the bilayer honey- 
comb optical lattice. In the low energy approximation, the dispersion relation has parabolic form 
and the quasiparticles are chiral. In the presence of the effective magnetic field, which is created 
for the system with optical means, the energy spectrum shows an unconventional Landau level 
structure. Furthermore, the experimental detection of the spectrum is proposed with the Bragg 
scattering techniques. 
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1 Introduction 

Recently, the studies of cold atoms in optical lattices are extensively developed. Optical lattices are 
crystals made of light periodic potentials that confine ultracold atoms[Tl [5]. Because of their precise 
control over the system parameters and defect-free properties, ultracold atoms in optical lattices 
provide an ideal platform to study many interesting physics in condensed matters |3| and even high 
energy physics 0]. 

Very recently, a strong interest has been raised in the two-dimensional honeycomb lattice f5l [6l [7] , 
for its physics is closely related to that of the graphene material [51 [HI UHl [HI [HI [ISJ, which has 
surprisingly rich collective behaviors. With the tight-binding approximation, graphene has a linear 
dispersion relation resembling the Dirac spectrum for massless fermions. In the presence of a magnetic 
field, it has the Landau energy level with square-root dependence on the quantum number n, instead 
of the usual linear dependence. In particular, the zero-energy Landau level exists at n = 0, which is a 
direct result of chirality. Recently, McCann et al. have studied the electronic states and unconventional 
Landau levels of the bilayer graphene arranged according to Bernal stacking fHl I15j . 

In this paper, we investigate the eigenstates and spectrum of ultracold fermions in the bilayer 
honeycomb optical lattice with a different stacking order from that in Reference [T31[TS] . In the absence 
of an effective magnetic field, the dispersion relation has parabolic form and the quasiparticles are still 
chiral like that in the monolayer system. In the presence of an effective magnetic field, which can be 
built by coupling the internal states(spin) of atoms to spatially varying laser beams [TOl [T71 [T51 [TOl . 
the spectrum shows an unconventional Landau level structure. The experimental detection of the 
spectrum is proposed with the Bragg scattering techniques. 
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2 The model 



We consider a system of ultracold fcrmions confined in the bilayer honeycomb lattice. The honeycomb 
lattice consists of two sublattices denoted by A and B. Then, the bilayer honeycomb lattice considered 
in our work is formed by coupling the B sublattices of the two layers with tunneling and leaving the A 
sublattices of the two layers uncoupled. One can create the bilayer honeycomb lattice in the following 
steps. First, one builds the monolayer honeycomb lattice as shown in Fig[T] (a) with three laser 
beams in the x — y plane and two laser beams along the z direction When the potential barrier 
of the optical lattice along the z direction is high enough, the vertical tunneling between different 
planes is suppressed seriously, then every layer is an independent two-dimensional honeycomb lattice. 
Secondly, one makes the triangular lattice as shown in FigH] (b) with red-detuning laser fields[22]. 
Finally, to realize the bilayer honeycomb lattice, one can put the triangular lattice and the honeycomb 
lattice together as shown in FiglT] (c) . There exists an additional micro-trap between sublattices B 
of every two layers in the honeycomb lattice. The additional micro-trap lowers the barrier between 
sublattices B of these two layers in the honeycomb lattice, or links sublattices B of these two layers 
in the honeycomb lattice as an intermediate point, so that sublattices B of these two layers in the 
honeycomb lattice are coupled. Following this scheme, many independent bilayer honeycomb lattices 
can be achieved (see FigH] (c)), so we only need to investigated one of them. For convenience, we 
assume that the whole system is trapped in a two-dimensional box, which can be achieved by adding 
four blue detuning endcap beams at the edges of the optical lattice in x — y plane 23J. With this box 
trap, the system can be considered to have the hard wall boundary condition approximately, so we 
can neglect the boundary effect in our discussion. 

In this scheme, the ultracold atoms have a A-type three-level configuration, the states |1) and |2) 
are degenerate states, which are assumed to be different Zeeman states on the same hyperfine level, 
and |3) is an excited state. The ground state |j) with j = 1,2 and the excited state |3) are coupled 
though two laser field with the corresponding Rabi frequencies VljC^'^^ , respectively [18]. The schematic 



representation of this scheme is as shown in FIG. [H The total Hamiltonian reads, H 
non-perturbative Hamiltonian is given by, 



Ho^Y.j dV^iiv) [-ri2vV2m + y(r)] ^„(r). 
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and the light-atom interaction Hamiltonian Hi is given by, 

Hi^ J dr 17ie*'^W'l'55(r)^i(r) + r22^'i(r)*2(r) +H.C. 
Diagonalizing the interaction Hamiltonian with the unitary transformation S, 
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The new field operators corresponding to the eigenstates are related with the old field operators as 
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In the new bases, and under the adiabatic condition ($1 jiJol'&j) ^ \Ei — Ej \ for j = 2, 3, we can apply 
the adiabatic condition and then neglect the populations of the states li'2) and l^a). Therefore, the 
effective Hamiltonian can be rewritten in the dark-state basis |$i) |16[ fTTl [TSl fT9l 1201 . 
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where A = — ?isin^ 9\/tp and 7i = 2^(— *?iV — A)^ + V"e//(r) is the single particle Hamiltonian with 
Ve//(r) being the effective trap potential. Here A is the effective gauge potential associated with 
the artificial magnetic field B = V x A. In the practical case, we choose two counter-propagating 
Gaussian laser beams as fl^e*'''' = fio exp[— (a; — a;j)^/(TQ] exp{~ikjy) {j = 1, 2), where the propagating 
wave vectors ki = —k2 = A:o/2 and the center position xi = —X2 = Ax/2 p^. Then the effective trap 
potential is, [18] 

, , fi^kl (1 + 1/4^2 fc2) 



and the effective vector gauge potential is 



fikt 



with d = (70/(4 Ax). Straightforwardly, one can obtain the effective magnetic field, [18] 

B = T Gz. (8) 

4dcosh2(x/2d) ^ ^ 

Practically, one may set d ~ 1mm and —O.Olmm < x < 0.01mm, so the condition \x/d\ <C 1 is 
satisfied. In this practical condition, the effective trap potential can approximately be written as 

which has an additional constant term compared with the original external trap potential. This 
additional constant term does not change the geometrical structure of the original trap potential, so we 
can drop out it as a constant chemical potential term. The effective magnetic field can approximately 
be written as 

B « B(o) + B(2) - - (10) 

where the quadratic term can be neglected for \B^'^^/B\ < 1.25 x 10"^ Thus, the effective magnetic 
field can be regarded as a homogeneous one in the regime considered in our scheme. For the typical 
parameter value ko ^ 2 x 10^m~^ and d ~ 10~^m, we obtain the magnitude of the effective magnetic 
field, B ~ 3.3 X 10-25J . s . m.-^. 

3 The effective low energy Hamiltonian 

Taking the tight-binding limit, we can superpose the Bloch states to get two sets of Wannier functions 
w^{r — Ti) and w^(r — r^) with a = 1,2, which correspond to sublattices A and B of layer a, 
respectively. In the presence of the effective gauge field we can expand the field operator in the lowest 
band Wannier functions as. 



*i(r) - E 
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Substituting the above expression into Eq.([5|), we can rewrite the Hamiltonian as H = Hq + Hi with 

a rieAj=l,2,3 
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and 



Hi = -t± ^[6l(r,)62(rO+H.c.], (13) 

where t is the tunneling parameter with t = — J drw^(r— ri)7ioWB(r— r^) ~— J drw*g{r—rj)'HoWA{j^— 
Ti); t± = — J drw*g^_^{r—ri)'Ho'WBi{r—rj); the energy shifts for sublattice A and B arc = Jo?rw^(r — 

ri)'HoWAir - r^) and es = J drwg{r - Vj)'HoWB{Y - rj) respectively, with Hq = -^V^ + Ve//(r). 
Here, for convenience, we can have dropped out a constant term in Hamiltonian (|12p . The three vector 
Sj in Eq. fT^ are Si = (0, — l)a, S2 = (\/3/2, 1/2) a, and S3 = (— \/3/2, 1/2) a, where a is the lattice 
spacing. The three vector Si{i = 1,2, 3) connect any site of sublattice A to its nearest neighbor sites 
belonging to sublattice B in every layer. 

Here, we assume the condition t± t being satisfied, so that we can consider Eq. ([13]) as a 
perturbation. We take the Fourier transformation to a(r) and b{r) as. 



a„(k) = e-'^^^^'aalr), (14) 

rieA 

6„(k) = ^^'•'■'"'^^(rO, (15) 

where r is the coordinate on x — y plane. Substituting the above expressions into Eqs. (fT2|) and ^3]) . 
we obtain the following Hamiltonian, 

^0 = 5]$]K(k)al(k)6„(k) +r(k)6l(k)a„(k)], (16) 

a k 

Hi = -t^ ^[6i(k)t62(fc) + S2(fc)tSi(fc)], (17) 

k 

where ^(k) is the single-particle energy spectrum without interlayer tunneling and defined via ^(k) = 

^ _2fk-s — — f^-' A.-dr) 

—tj^j^i 2 3'^ ' Jo ', The energy spectrum contains two zero-energy points at K± = 
±(47r/3-\/3a, 0) around which it is linearized. Neglecting the coupling between the Fermi points K-t, 
the total Hamiltonian H can be expand around the contact point K_|_(K_) in coordinate space. 
Without loss of generality, we expand the total Hamiltonian around the contact point K+ as[T3[T3], 

H = [ d^rij\v)'Hij{v), (18) 



where the spinor -0 — (-01 "01 "02 riY' f^'' the Dirac point K+. Here, 7i takes the 4x4 matrix form, 

(19) 



/ vpn^ \ 

wftt —t± 

vpn'^ 

\ -tj_ VpTT / 



where tt = -frx + iTTy and tt^ — fcx with tTx — Px Ax/fi and T^y ^ Py — Ay/h, and vp = 3at/2h 

is the Fermi velocity. Here, t± t is assumed. Eliminating the dimer state components ifji and 
we can reach a two-component Hamiltonian describing effective hopping between the A1-A2 sites 

n..^f( I ^'M, (20) 

2m \ TT'n J 

where m = t±/2vp. 
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4 Energy bands 



First, we consider the case without gauge fields, i.e. -kx — —idx and r^y = —idy. The eigenfunctions 
of the Hamihonian ([^0]) are given by 

/,k(r) = ^exp(zk.r)(^ J (21) 

where is the area of the system, and s denotes the conduction band with s = +1 and the valence 
band with s = — 1. The corresponding eigenenergies are 



where k = ^k^ + ky. This dispersion relation has parabolic form as shown in FiglS] Here, the 

quasiparticles are still chiral like that in the monolayer system. The pseudospin vector n = (1, 0) is a 
constant for any wave vector k in our work, while n — (cos(20), sin(20)) for k = (fc cos (p, k sin cj)) in the 
bilayer graphcnc with Bcrnal stacking order as in References |141I15] . Thus, in our bilayer honeycomb 
lattice configuration, the Berry phase 27r in the bilayer graphene with Bernal stacking order is absent. 



5 Unconventional Landau levels 



For the case with an effective magnetic field in the Landau gauge (0, Bx, 0), the eigenfunctions of the 
Hamiltonian (1201) can be obtained as, 
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with sgn(rt) = (1, 0, —1) for (n > 0, n = 0, 7i < 0) respectively, for n ^ 0, and 
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for n = 0. Here, 



are harmonic oscillator eigenstates as 
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where the quantum nubmer n is an integer and u = ^fijB. The corresponding Landau energy levels 
are En = nhujc with n = • • • , —2, —1, 0, 1, 2, • • • and ujc — 2vpB/t±, which is shown in Fig[4l This 
spectrum has a linear form like the conventional Landau level spectrum. However, the quasiparticles 
are chiral in our scheme. The zero modes exist and the zero energy level is twofold degenerate 
compared to the non-zero energy levels. 

To give a numerical evaluation, the typical values of the parameters can be taken as t ~ 10~^°J, 
tj_ ^ a ^ 200nm, B ~ 3.3 x 10~^^J • s • m~^. We can then estimate the magnitude of the 

cyclotron frequency ujc — 5.7 x lO'^s^^ and obtain the first gap of Landau level for the monolayer 
honeycomb lattice A ~ 6 x 10^'^^ J. The temperature required to keep atoms in the zeroth Landau 
level is 43 nK. 



6 Bragg spectroscopy 

It is not easy to measure the Hall conductivity of cold fermionic atoms in the bilayer honeycomb lattice. 
However, an available method to detect the unconventional Landau levels of ultracold fermions on 
the bilayer honeycomb lattice is the Bragg spectroscopy [25 , which is extensively used to probe the 
excitation spectrum in condensed matter physics. In the Bragg scattering, the atomic gas is exposed to 
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two laser beams, with wavevectors ki and k2 and a frequency difference w. The light-atom interaction 
Hamiltonian can be written as, 

Hb= f^i3e'i"-|/,„k+q)(/.,.k|+H.c., (26) 

Sl,S2,k,q 

where Si = 1,2. In our case, we consider the case of half filling, i.e., the bands with n < are 
fully occupied and the bands with n > are empty. The half filling state can be prepared with the 
coherent filtering scheme proposed in Reference [5S]. From the Fermi's golden rule, we obtain the 
dynamic structure factor as follows, 

a 

x6{nuj - Ef} + Ea) (27) 

where N is the total number of atoms in the system; denotes the initial state and (pa^ denotes 
the final state; a represents all quantum parameters of quantum state. 

For simplicity, we assume that the direction of q is the same as the one of y axis in the wavevector 
space, i.e. q — qey. Following the above formulae, we can straightforwardly evaluate the dynamic 
structure factor ^(q, w). FigO^a) shows the dynamic structure factor S{q,uj) as a function of ui for 
the case without an effective magnetic field. We can find that S'(q, lu) is zero when lo is under 
and is finite constant for uj above h^q'^/2m. Figl5](b) shows the dynamic structure factor S{q,uj) as a 
function of uj with q = l.Q^J B /% for the bilayer honeycomb lattice. Similarly, when a zero-level atom 
is excited to other states, the peaks are obtained aiuj— (1, 2, • • • , n, • • ■)K^t^ /1it±^ with A = ia^J B /2fL^ 
which are marked with red stars in FIGE] (b). The distances between the neighbor peaks marked with 
red stars in Figl3] (b) are identical. 

7 Conclusion 

In summary, we have proposed a scheme to investigate ultracold fermionic atoms in the bilayer hon- 
eycomb lattice for the cases without and with an effective magnetic field. The effective magnetic field 
can be built with optical techniques. For the case without an effective magnetic field, the dispersion 
relation has parabolic form and the quasiparticles are chiral. For the case with an effective magnetic 
field, there exist unconventional Landau levels that include a zero-mode level. We have calculated the 
dynamic structure factors for the two cases and proposed to detect them with the Bragg spectroscopy. 
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Figure 1: (a) The independent monolayer honeycomb lattice, (b) The adding triangular lattice, (c) 
The bilayer honeycomb lattice built with putting (a) and (b) together. 




Figure 2: The light-atom interactions between fermionic atoms and two laser beams. 
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Figure 3: Energy bands of cold fermionic atoms in the bilayer honeycomb lattice without an effective 
magnetic field. 
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Figure 4: Landau levels of cold fermionic atoms in the bilayer honeycomb lattice with an effective 
magnetic field. 
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Figure 5: The dynamic structure factors S{q, lo) (not scaled) for cold atoms in the bilayer honeycomb 
lattice in the cases (a) without an effective magnetic field and (b) with an effective field. Here, uj is 
in units of 9at^q^/2ht± with q = l.Q^B/h. 
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